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Abstract 

Using a variational model based on non-linear elasticity we investigate whether in 
a cubic-to-tetragonal phase transformation it is energetically preferable to nucleate 
martensite within austenite. Under minimal growth assumptions on the free energy 
density W away from the wells, we derive explicit upper bounds on W qc (T), i.e. on the 
macroscopic free energy density of a region that has been macroscopically deformed by 
the identity map. The bounds only depend on material parameters, the temperature 
difference and the growth of W away from the wells. By comparing W qc ( I) and W(V) 
we conclude that nucleation is always energetically preferable and are able to provide 
quantitative upper bounds for the gain in negative energy due to nucleation. 
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1 Introduction 

The high growth velocity of heterogeneous temperature induced nucleation of martensite, 
which is of the order of magnitude of sound ( [BM53] . |Chr02l p.977 ff.]), is one of the 
reasons why the process of nucleation is not yet satisfactory understood. As pointed out 
in |Chr02l Ch. 23] “the kinetics of martensitic transformations are depended mainly on the 
nucleation process as each plate grows rapidly to its limiting size” and as a consequence “it 
is often found that the plates of martensite formed (...) have a constant lenghtThickness 
ratio” until “growth ceases because of plastic deformation”. Different models (e.g. [ BKS1L 
ICoo73L IOC791 IOR92j ) have been proposed to explain the formation and predict properties 
of such plates. It has been shown theoretically and experimentally that the availability 
of a free surface can have a significant influence on the nucleation process. In [KPG89j a 
single crystal of Fe-Ni was coated with a Ni-rich surface layer to inhibit nucleation from 
the free surface with the effect that the formation of martensite was completely suppressed 
during slow cooling. Applying a similar procedure to an In-Cd alloy however did not show 
a difference between the coated and uncoated crystal. In a related experiment on Cu- 
Al-Ni single- and polycrystals with varying grain-sizes it has been shown in |SC79| that 
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the energy needed for bulk nucleation is significantly higher than that required for corner 
nucleation or to move a planar interface. A single crystal of the same alloy was considered 
in [BKSll j where it was shown experimentally and theoretically that corner nucleation is 
energetically preferable over interior nucleation. However if a free surface is not available, 
for instance because of a relatively large grain-size or by a coating as described above, it 
is of interest to establish criteria for interior nucleation. 

The present paper investigates the energetic properties of interior nucleation using a 
variational model based on non-linear elasticity first introduced in |B.187| . In contrast to 
|OR92] . where a nucleus forms with an array of interfacial dislocation, in this model the 
material is assumed to deform purely elastically without considering plasticity. Further 
interfacial energy is not taken into account. This is justifiable by recent experiments in 
j DSjSll j on Ti-Ni-Pd and Ti-Ni-Au crystals with very low hysteresis where it has been 
observed that structures with very complex interfaces form, even though less complex 
interfaces are already kinematically compatible, so that the contribution from interfacial 
energy must be negligible. 

A particularly interesting case of interior nucleation is self-accommodation. Using the 
Ball-James ([ BJ87 ]) model it has been shown in (Bha92l that a necessary and sufficient 
condition for self-accommodation is that the austenite and martensite have equal volume. 
In this case the martensite can rearrange itself such that it matches the austenite along 
any surface and thus interior nucleation is energetically as preferable as nucleation from 
the boundary. However in general we expect to have a volume change during transforma¬ 
tion, so that any inclusion of martensite within austenite necessarily has a region where 
it is neither on the austenite nor on the martensite wells. We will refer to this region as 
the interpolation layer and it can loosely be understood as the elastic equivalent to having 
interfacial energy and plasticity. By the classical theory of heterogeneous nucleation in 
order to transform a given region of austenite to martensite “the negative free energy 
change resulting from the formation of a given volume of a more stable phase [the marten¬ 
site] is opposed by a positive free energy change due to the creation of an interface [the 
interpolation layer] between the initial phase and the new phase” f |Chr021 p. 4]). Thus our 
criterion for interior nucleation is that the energy of the inclusion is lower than the energy 
of the parent phase, i.e. the austenite. It seems to have first been shown qualitatively in 
(Ball3| that for a sufficiently flat cylindrical inclusion nucleation is always energetically 
preferable over the pure austenite. In this paper we develop quantitative estimates on the 
energy difference between a nucleus and the austenite depending on the temperature and 
the assumed growth of the free energy density away from the wells. 


Nomenclature 

W 1,00 (Q) 

Space of continuous functions from H to M 3 with bounded derivative. 

w 0 1,oo (fi) 

VF 1,00 (H) with zero boundary conditions. 

r 3x3 

Space of 3 x 3 matrices. 

SO(3) 

Space of rotation matrices, i.e. R e M 3x3 such that R T R - I, deti? = 1. 

p 24 

Crystallographic point group of a cube. 
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2 Description of the Model 

Following |BJ87| microstructures correspond to minimising sequences of the total free 
energy 

E (y) ■■= f n w e(Dy(x))dx, (E) 

where !lcK 3 represents the undistorted austenite at temperature 6 and y e VF 1,0 °(f2) 
maps a point x in the reference configuration to a point y(x) in the deformed con¬ 
figuration. In this continuum model the mesoscopic free energy density W depends 
only on the deformation gradient F = Dy(x) and temperature 9. The deformation 
y(x) = x , Dy(x) = I corresponds to pure austenite and thus, since we are interested in 
modelling inclusions of martensite within austenite, y has to satisfy the homogeneous 
displacement boundary condition y(x) = x on dQ. Furthermore W respects frame indif¬ 
ference, i.e. W(RA ) = W(A) for all R e 50(3), A e M 3x3 and the cubic symmetry of the 
austenite, i.e. W(AQ) = W(A) for all Q e V 2i . Let K m = SO(3)OilP 24 = Uf =i SO(3)I7i, 
U\ = diag(r/ 2 ; Vi), U 2 = diag(? 7 i, 772 , yi), U 3 = diag(r/i, 771 , 772 ) denote the set of martensitic 
wells. Then for a fixed temperature 6 * below the transition temperature we assume 

W(R) = 0 for R € S0(3), W(U) = -r for U e K m , W > -r on R 3 x 3 \K m . (1) 

Clearly the energy functional (IE]) is bounded from below by F(y) > -r Vol(P). Hence any 
sequence {yj} c H /1 ’ 00 (P) such that E(yj ) -»■ -tVoI(P) corresponds to a microstructure 
that can only contain variants of martensite and thus can only exist if the volume does not 
change during transformation. If the transformation does involve a change in volume, i.e. 
det U\ ± 1, it is of interest to find the minimal value of the free energy E. Since the energy 
is proportional to the volume of the inclusion the relevant quantity is the macroscopic free 
energy density, i.e. the average energy per volume element, denoted by W qc . W qc (l) is 
the free energy density of a region that has been macroscopically deformed by the identity 
map, i.e. that obeys the boundary condition y(x) = x. A characterisation (see |Dac82| l of 
W qc is given by 

W<1C W = 7717m inf [ n W(I + D<Kx))dx, 

Vol(H) <fxWQ } °° (f 2 ) ~>n 

so that inf E(y) = Vol(P)VF 9C (lI). In particular whenever W qc (I) < 0 there exists a defor¬ 
mation y(x) = x + 4 >(x) that is energetically preferable over the pure austenite. 

3 Interior Nucleation of Martensite 

In the following we will always assume det U\ + 1 so that we cannot have stress-free inclu¬ 
sions. The goal is to construct a test function <p such that the corresponding deformation 
y(x) = x + 4 >(x) has negative total free energy and is thus energetically preferable over the 
pure austenite corresponding to y(x) = x. We will explore two different approaches. In the 
first approach we only allow deviations from the austenite and we show that the geometry 
of the inclusion plays a crucial role. In the second approach we only allow deviations from 
the martensite wells and use the principle of self-accommodation to derive our estimates. 

Deviations from the Austenite 

By [BJ87j it is known that a fine mixture of two rank-one connected martensitic variants 
can form a stress-free planar interface with the austenite. That is, there exists a rank-one 
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matrix b ® m, |m| = 1 such that 

XRUi + (1-X)RRU 2 = I + b®m, (2) 

where RU\ , RRU 2 e K m such that U\ - RU 2 = a® m, |n| = 1 and A e (0,1). We consider 
an inclusion with a core of pure martensite. Since the volumes do not agree we have a 
transition layer where the deformation gradients cannot be supported on the wells. We 
show that the energy penalty from the transition layer can be overcome by the gain in 
negative energy due to the presence of the martensite as long as the inclusion is flat and 
thin. 

Definition 1. (Cylindrically symmetric inclusion with a simple laminate core) 

Let n = n 1 un 2 be such that the core fli consists of a simple laminate of martensite that 
is compatible with the austenite in the sense of ([ 2 ]) . By rotating the entire system we may 
assume without loss of generality m = e 3 . We define the two regions as follows 

:= {x e M 3 |x 3 e [0, ip(r )]}, fl 2 ■= {x e M 3 |x 3 e (tp(r), cr^(r)]}, 

where r 2 = xf + x|, a > 1 and 0 < ip(r) e LL 1,oo ([0, 00)), lirn^oo ip(r) = 0. Note that we 



assume symmetry with respect to the x\ - x 2 plane and only define the upper half of the 
inclusion. We further assume that Vol(fli) = 27 t / 0 °° r<p(r)dr = 1 and thus Vol(ri 2 ) = o - 1. 
We define our test function <f> e lL / 0 ’ DO (n) in the regions Hi and Ll 2 separately. In Hi we 
set 4>i(x) = </>[ n := b ■ x 3 and in Ll 2 we interpolate linearly in x 3 direction from b ■ <p(r) to 
0 by defining 

Mx) = (t>\n 2 : = ~ ( r )) ’ 

so that (j) = (j)\ + (j ) 2 € W 0 1 ’°°(n). 

Theorem 1. Let W be as in (jT]) and assume the following upper bound on W away from 
the austenite 

W(I + H) <C\H\ 2 . (3) 

Then the inclusion given by Definition^ provides the following upper bound 

W qc (T) < -r (l + 2/3 + 2^/3 + /3 2 ) _1 , (4) 

where fd ■= ^|b| 2 . 
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Proof. In Hi the gradients of f> are supported on the martensite wells and therefore E{(f> i) = 
-rVol(Hi) = -r. In Q '2 we estimate the energy using the upper bound on W {I + H). 
We have Df> 2 (x) = g Vri(x) with Vu(x) = l) . Noting that 

Df> 2 (re) is independent of X 3 and using the coarea formula (see [Fed96l Ch. 3]) we arrive 
at 


E(<t> 2 ) < 


C 


(a - !) 


Ib| 2 + 


a 


(u- !) 


7, 


where 7 := C 7 r|b | 2 rip(r)(p' (r) 2 dr. Combining the two estimates we obtain 


E(4>) < -t + 


c K- " 


(cr-1) (cr - 1) 
and thus for 7 -> 0 and since Vol(H) = a we have 


7 


W qC {l)<T\--+ C 1 


a t a (a - 1 ) 


|bf 


1 


r — + 


P 


a - 1) 


where /3 ■= ^|b| 2 . Minimisation with respect to a gives the desired bound. □ 

Remark 1. We note that we derived this upper bound by taking the limit 7 -> 0 and 
thus /g 00 r(p(r)(p'(r) 2 dr -*■ 0 which corresponds to <p'(r) -> 0 as the volume is fixed. Thus 
within the class of inclusions under consideration the lowest energy is achieved by a very 
flat and thin inclusion. 


We can derive the same upper bound on W qc { I) by estimating W rc ( I), where W rc 
denotes the rank-one convex envelope of W (see e.g. I)ac()7 ). The proof uses the general 
result that W qc < W rc and is simpler than the proof of Theorem [I] but gives no information 
on the structure of the inclusion. 

Lemma 1. Under the same assumptions as in Theorem^ 7] an upper bound on W rc (ll) and 
thus on W qc (V) is given by 

W rc {l) < -r (l + 2/3 + 2y//3 + /3 2 )~ i , 


where j3 ■= y|b| 2 . 

Proof. An upper bound on W rc (see |Dac071 Thm. 6.10]) is in its simplest form given by 
W rc (A) < inf{AW(Ai) + (1 - A )W(A 2 ) ■ XAi + (1 - A )A 2 = A,rank(Ai - A 2 ) < 1}. (5) 


The estimate can be refined by applying the same estimate to each of the terms W(Af), i = 
1,2 on the RHS. We now decompose I = A(I + b <g> m) + (1 - A)(I + eb g> m), £ = ^ and 
II + b ® m = XR 1 U 1 + (1 - X)R 2 U 2 , Ri, R 2 e SO(3) such that rank(Ri!/i - R 2 U 2 ) = 1. Then 
by using ([ 3 ]) and ([5]) we have 


W rc (l) < XW rc (l + b 0 m) + (1 - A)W(I + b 0 m) < r 

A - 1 




with fj := ^|b| 2 . Minimisation with respect to A gives the desired estimate. 


□ 
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Deviations from the Martensite 


In many materials detUi « 1, so that matrices in a relatively small neighbourhood of the 
martensite wells have unit determinant. Allowing deviations from the martensite wells we 
now consider an inclusion with a core of variants that are close to the martensite wells 
and allow self-accommodation and thus do not require a transition layer. 

Theorem 2. Let |det t/ -1 / 3 - 1| « r/i - 1|, \r ]2 - 1|, U e K m = Uf =1 SO(3 )Ui and let W be 
as 0 with 

W{V) < W(U) + C dist 2 (V, K m ) (6) 

for all V e M 3 * 3 such that dist(V, K m ) < |1 - det U~ x ! 3 \. Then 

W qc (I) < W rc (I) < -r + C(det U- 1/3 - 1) 2 \U\ 2 . (7) 


Proof. By Bhattacharya’s construction (see jBh a92j ) in a cubic-to-tetragonal phase trans¬ 
formation there exist laminates with macroscopic deformation gradient I iff dett/i = 1. 
Therefore let us set Ui := (det U)~ l ! 3 Ui so that det Ui = 1. By 0 we have 

W(Ui) <-T + c(l- det tr 1/3 )" \U \ 2 

and thus using 0 we can estimate 

W qc (I) < W rc ( I) < W(U) + C(l- det U- 113 )' \U\ 2 . 


:□ 


Remark 2. We note that unlike Theorem [l] this construction is only preferable over the 
pure austenite for r sufficiently large, corresponding to a temperature 9* sufficiently far 
below the transition temperature 9 C , or for det U sufficiently close to 1. A comparison with 
the estimates in Theorem [I] reveals that ([l]) provides a better bound than Q for small 
volume changes and a worse bound on W qc (I) for large volume changes. For det U -*■ 1 
we recover the case of self-accommodation, i.e. 

-T < W qc (I) < -T + lim C(detcr 1/3 - 1) 2 |?7| 2 = -r. 

det C/-»l 

Remark 3. For simplicity the present paper only considered a cubic-to-tetragonal phase 
transformation. However, in the first part the only assumption on the material was the 
existence of an interface between a microstructure of martensite and austenite. In the 
second part, the only assumption on the material was that it almost satisfies a criterion 
of self-accommodation (e.g. in a cubic-to-* transformation this holds iff det Ui ~ 1, cf. 
[Bh a92] ). Thus as long as either or both assumptions are satisfied all respective conclusions 
remain the same. 


4 Conclusions 

Without making any a priori assumptions, apart from cylindrical symmetry, on the shape 
of an inclusion of martensite within austenite we have shown that - within the class of 
inclusions under consideration - the lowest energy is achieved for a flat and thin shape 
as observed in experiments. We showed that W qc (I) is always negative, so that within 
the model, it is always energetically preferable to nucleate martensite within austenite. 
Further by comparing two different types of nucleation in Theorem [l] and Theorem [2] we 
could show that for materials that almost fulfil the criteria for self-accommodation, i.e. 
detC/i » 1, it is energetically preferable to form a self-accommodating inclusion rather 
than an inclusion with a core of a simple laminate. 
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